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This paper is concerned with the out-of-equilibrium two-lead Kondo model, considered as a model 
of a quantum dot in the Kondo regime. We revisit the perturbative expansion of the dot's magneti- 
zation, and conclude that, even at order in the Kondo interactions, the magnetization is not given 
by the usual equilibrium result. We use the Schwinger-Keldysh method to derive a Dyson equation 
describing the steady state induced by the voltage between the two leads, and thus present the 
correct procedure for calculating perturbative expansions of steady-state properties of the system. 
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I. INTRODUCTION 



In recent years, much experimental and theoretical work has been devoted to exploring the properties of so-called 
'quantum dots' M. These are mesoscopic devices in which a 'dot' containing a small number of electrons is isolated 
from two macroscopic leads (denoted 'left' and 'right') by potential barriers, through which tunneling processes take 
place. Experimentally, these systems are small devices fabricated using a two-dimensional electron gas ||, |J, [f| 
or carbon nanotubes ||. The number of electrons on the dot, N, is controlled by a gate voltage V g . From the 
■ experimental point of view, one is primarily interested in the current flowing through the dot as a function of V g 
and of the potential difference V between the two leads. When tunneling is weak enough, the Coulomb blockade 
phenomenon appears Q : the conductance through the dot is essentially zero except in the vicinity of certain special 
values of V g , where the energy difference between the ground states of the dot with N and N + 1 electrons vanishes. 
£C) , At these points, conductance peaks are observed. 

A simple model for such a system is the Anderson model, where the localized level represents the dot and the 
hopping term describes its hybridization with the leads. In this paper, we shall restrict ourselves to the regime where 
the occupation of the dot is not fluctuating, and where naive application of Coulomb blockade ideas would predict a 
strongly suppressed conductance. However, when N is odd there is exactly one unpaired spin, which is coupled to the 
leads via a Kondo interaction. In this case, in the linear response regime, spin physics opens up a new transmission 
' channel via the Kondo effect. Although the occupancy of the dot remains fixed, spin-flip interactions permit the 
formation of strong dot-lead hybridization for temperatures T < Tk, where Tk is the Kondo temperature. As the 
temperature approaches zero, this leads to unitary limit conductance (G — Go = 2e 2 /h) via a sharp resonance at 
& \ the Fermi surface — the Abrikosov-Suhl resonance. This effect was predicted in the context of quantum dots fourteen 
years ago [0, || , and was recently observed in a series of experiments ||, ||, ||, [|, ||] . 

By contrast, in the large- voltage regime, full non-equilibrium calculations are required, and much less is known than 
close to equilibrium. The out-of-equilibrium Anderson and Kondo models have been studied by various methods in 
the last decade. Much of this work j^, [II], [llj has concentrated on the non-crossing approximation (NCA) approach, 
£j , which adopts a slave-boson description of the problem, and then renders it tractable by neglecting certain vertex 
corrections. This yields a picture in which the Kondo resonance in the density of states is both split and broadened 



as the voltage V between the leads is increased. Recent work 12 has explored a new approach, where one attempts 
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to use the Bethe Ansatz results for the Anderson model to construct a Landauer-type picture of transport through 
the dot. This approach also involves approximation, when one comes to construct the 'in' and 'out' scattering states 
from the dressed excitations of the model. Another thread [[l3| 14 has involved studying the Anderson model via 
perturbation theory in the on-site Coulomb repulsion U . While these works provide approximate information on the 
behavior of the current- voltage characteristic, they shed little light on the nature of the many-body state of the system 
when V > Tk- 

In particular, a basic question recently debated |0, |l5| is whether the Kondo problem has a strong coupling regime 
at low temperature and high voltage. In discussing this point, a previous paper p5[ used a second order perturbative 
expression for the magnetic susceptibility in the out-of-equilibrium steady state induced by V. This putative result 
was however incorrect, even at order in J, the strength of the Kondo couplings; in this note, we correct this result at 
order 0, discuss the related physics and present a systematic method for calculating higher-order corrections in J. Our 
main conclusion is that, even at order in J, the Kcldysh function of the spin deviates from its equilibrium value and 
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therefore the steady state magnetization of the dot is not given by M eq = ^ tanh(£?/2T). Rather, it must be computed 
by solving a transport equation. We emphasize that this issue is not directly related to the so called "decoherence 
time" but is a basic point about perturbation theory to be addressed before discussing out-of-equilibrium RG 
equations and the existence of a strong coupling regime at large voltage. In particular, it has ramifications f or othe r 
physical quantities, e.g. the current and the current-current correlation function in a magnetic field (see section [II C 2] ). 



The paper is organized as follows. In Section II, we present the model, our main result (the magnetization at order 
given by Eq. ||), and the associated physical discussion. In Section III, we present a detailed pedagogical derivation 



of (^) using the Keldysh method. Finally, in Appendix [d[ we present more details about the computation presented 
in ref. p"5| , and explain why it was incorrect. 

II. RESULTS AND DISCUSSION 

Our starting point is the Hamiltonian of the two-lead Kondo model. For a discussion of the modeling of the 
quantum dot, and for a derivation of this model from the Anderson model via an out-of-equilibrium Schrieffer- Wolff 
transformation, we refer the reader to the literature, in particular pl| and references therein. The Hamiltonian is: 

H =^2(e k - ^ a )c f aka c aka + H IeB + H tlans - B S S Z - B c crc^ k(T c ak(T , (I) 

aker ak<r 

Hrcfi =Jr ( C flk<T^' C flk'^') ' $ + 
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where c^ kcr creates an electron in lead a € {L, R} with momentum k and spin a, and Jl, Jr and Jlr = (Jrl)* are 

Kondo coupling constants between the electrons and the spin of the dot S. The first term in H describes the electrons 
in the leads, with £k being the bare energy of an electron of momentum k at zero voltage (the same for each lead) and 
fi a is the potential in the lead a. Each lead consists of a free electron gas with a density of states p(e) of bandwidth 
D: ultimately, we will be interested in the result in the large bandwidth (D — ► oo) limit, but the computations are 
first performed for finite D. We will make the physical assumption that the leads are in thermal equilibrium at a 
temperature T. The voltage is applied by taking the chemical potentials of the two leads to be different, — Hr = V . 
-ffrefl describes regular Kondo processes, where an electron from a given lead is spin-flip scattered back into the same 
lead; -fftrans describes 'spin-flip cotunneling', where an electron from one lead is spin-flip scattered into the other lead. 
If the model (j|) is derived from the Anderson model, one finds that JrJl = \Jrl\ 2 - In this paper, however, we 
relax this relation between the coupling constants and treat them as independent parameters. Finally, the last terms 
represent the coupling to the magnetic field. We allow two different magnetic fields, B s for the spin and B c for the 
lead electrons; this permits the calculation of the total and local spin susceptibilities within the same computation. 

In this paper, we shall be interested in the values taken by the following quantities in the non-equilibrium steady state 
induced by the voltage V: the dot magnetization Md ot , the magnetization of the leads Mi oa d s , the total magnetization 
Mtot, and the total and local susceptibilities. These are given by: 

M dot (B s ,B c ) = {S z ), (2a) 

Meads (-Bs, -Be) = ( CTC iko- C " kCT 

(2b) 

\aker / 

M tot (B) = M dot (B, B) + MieadsCB, B) - A/ Pauli , (2c) 
dM tot (B 



Xtot 



Xloc = 



dB 

dM dot (B s ,B c = Q) 



(2d) 



dB s 



(2e) 

B,=0 



where angle-brackets (...) denote an expectation value taken in the steady (i.e. long-time) state of the system. Mp au u 
is simply the Pauli paramagnetic contribution from the lead electrons which would be present even in the absence of 
the impurity, and which we therefore exclude from M to t- We consider the perturbative expansions of these steady 
state quantities; more precisely, we define: 

Jr = 9 r J J l = 6lJ Jrl = 0rlJ (3) 
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and we let J go to zero while keeping the coefficients Or, Ol and Orl fixed. In the following, the expression "order 
n" refers to the order n of this expansion in J. 

Our main result is the order term of the perturbative expansion of the magnetization: 

1 ,/b\ HI)! 1 



2 V 2T/ l 

2 



M tot (B) = M dot (B, B c ) + O(J) = £ tanh ( ^ ) -= 1— ^ ™ RL / 2 ^ + 0( J), (4) 



where </? is defined by: 



(At this order, Mdot{B, B c ) does not depend on £? c .) As a result, the magnetic susceptibility at order is: 



and in particular for |V| — * oo 

1 + 



X(T,V) = ± - — ■+()(■/.. ,(,) 



A striking feature of (Q) is that the magnetization at order in J is not given by the equilibrium expression 
M cq = itanh( J B/2T). This may seem surprising: since the couplings to the leads are relaxed to zero, why should we 
not find M cq , the magnetization of a free spin? Physically, the answer is that M cq is not the "magnetization of a free 
spin", but rather the magnetization of a spin weakly coupled to an equilibrium bath. On the other hand, Eq. gives 
the magnetization of a spin weakly coupled, via the Hamiltonian (|l|), to two leads with different chemical potentials. 
It is simply an out-of-equilibrium extension of the Curie law. In particular, we emphasize that the finite susceptibility 
at T = cannot be interpreted as a renormalization (or a screening) of the spin due to the voltage. 

At long times, the state of the spin is completely determined by the properties of these leads, and hence so is 
its distribution function, which describes the population of its two states as a function of temperature, voltage, 
and magnetic field. Since the whole system is not in equilibrium, this steady state is not described by the Gibbs 
distribution; in particular, the fluctuation-dissipation theorem (FDT) need not hold. Hence the magnetization need 
not be (and is not) M eq at order in J. Rather it should be computed by solving a transport equation in the 
steady state regime, i.e. a quantum Boltzmann equation. (At dominant order this procedure is equivalent to using a 
semi-classical master equation; see Appendix |a|.) 

The crucial point is that Jr, Jl, Jrl are relaxed to zero, assuming that they are still bigger than the coupling of 
the spin to any other thermal bath. If we were to take into account such a coupling (denoted by Jo) then the result 
would cross over to M eq when Jr, Jl, Jrl <C Jo (all the couplings going to zero while maintaining fixed ratios). In 
fact, (|j) implies that the equilibrium value is only an upper bound: < Af to t(J = 0)/M eq < 1, which follows from 
the convexity of ip. This bound is saturated only in equilibrium, i.e. for 8rl — * or V — ► 0. 

Moreover, as expected on physical grounds, the non-equilibrium result is much less universal than the equilibrium 
one. In equilibrium, the magnetization of a spin weakly coupled to a bath depends neither on the properties of the 
bath except the temperature T nor on the form of the couplings between the spin and the bath. Neither of these 
statements holds true for the non-equilibrium Kondo model: Afdot (at order 0) depends not only on T but also on the 
voltage V and on a ratio involving the parameters Or, 9l, and Ojil. Thus M do t is perturbative in J, but not in the 
three variables Jr, Jl, Jrl, and in particular it is not analytic in Jr, Jl, Jrl around (0, 0, 0). One can find such a 
dependence on the ratio of couplings even in a simple free model (a single level coupled to two leads), as illustrated in 
Appendix [b|. Furthermore, it should be remembered that we computed AJjot with free leads: introducing interactions 
in the leads would change the function (p (the important quantity being the electron-electron bubble in the leads). 
For example, we expect a different result to hold for a spin coupled to Luttinger liquids, even at order 0. 

It should be noted that the result (^) gives a non-trivial expression for the magnetization even at zero temperature. 
Since 4>(x) — > |x| as x — > oo, we find that 

1 I S R + S L 
1 1 ^32 



B 1 + 

2 [i(\B + V\ + \B-V\) + ^\B\ 



M tot (B, V)\ T=Q = — I — - - SL. ———— | + (J). (8) 
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We then have two cases (we can take B,V > 0): for B > V > 0, we obtain the equilibrium result M to t(B,V) 
sgn (B)/2 + O(J), but for < B < V, the magnetization is a still a function of B/V: 




Aftot L n » — — *± , (9) 



in agreement with (Q) in the limit \V\ 3> \B\. Physically, there are two sources of energy available to flip the spin: the 
thermal fluctuations of both baths (represented by T) and the fact that an electron can go from L to R and give to 
the spin an energy of order V. If we decrease B from high values (B 3> T,V), the spin is locked until B reaches the 
largest of these energy scales. Thus the magnetization at zero temperature is expected to saturate only for B > V. 
Similarly, the susceptibility is in general expected to behave as x ~ ^/B, where E is the largest energy available to 
flip the spin. The fact that our result is still non-trivial at zero temperature implies that it could be seen in numerical 
computations, such as the DMRG approach of Cazalilla and Marston p7| . 

III. DERIVATION OF EQ. 

The purpose of this section is to present the derivation of Eq. (^) and more generally the procedure for obt aining 



the perturbative expansion of physical quantities in the steady state. It is organized as f ollow s: in para graph III A 



we first give an overview of the derivation; full details are given in the following sections ( III B and QIC) and in the 
appendices, including a presentation of the Keldysh method. 

A. Overview 

When doing perturbation theory in the steady state, there are two important small couplings: J, the strength of 
the Kondo couplin gs, and s, the small regulator that appears in Green's functions such as the "spin" retarded function 



(defined in section III B 3) 



R%{oj) = i— . (10) 

u) — B + is 

The scale s should be thought of as being due to the coupling to an auxiliary thermal bath, which in the physical 
system would be the substrate. If the impurity is coupled to such a bath whose density of states is by a coupling 
g, one finds s ~ g 2 P2 (in the large bandwidth limit). In the physical quantum dot system, J»s, meaning that the 
correct order of limits to take is s — > followed by J — > 0, as pointed out in section II. 

However, 'straightforward' perturbation theory in J takes the limit in the opposite order: one first expands in J 
while keeping s finite, and only then takes s — > term by term in the perturbation series. In equilib rium, these 



two limits commute, but out of equilibrium they do not. This is explicitly shown in section [II CI; the out-of- 
cquilibrium Keldysh Green's function is not analytic around (s, J) = (0,0). This non-analytic behaviour is quite 
generic in situations where an impurity is coupled to several leads, and may be seen even in a simple free model (see 
Appendix |b|). The signature of that non-commutativity is that 'straightforward' perturbation theory fails: its terms 
(starting at 0(J 2 )) exhibit divergences of a 1/s form. These are similar to the infrared divergences in equilibrium 
perturbation theory that signal an incorrect choice of reference state. As shown below, 1/s divergences at order J 2 
signal an incorrect choice for the Keldysh Green's function at order 0, or equivalently of the distribution function 
which describes the non-thermal population of the two levels of the spin. Consequently, contrary to the claims of 
p5| , there is no possibility of regulating these 1/s divergences order by order in J. We discuss the putative regulation 
procedure of |lq ] in Appendix and explain why it is incorrect. 

The solution is to begin with the Dyson equation in the steady state, considered as a functional equation for the 
full Green's function G, using the skeleton self-energy diagrams (see Appendix [I] for further details). After taking 
the s — > limit in this equation, a perturbative expansion for G may be inserted into it, and a solution obtained 
order by order in J. Solving the Keldysh component of this Dyson equation is equivalent to solving the quantum 
Boltzmann equation; as shown below, the solution is non-thermal even at order in J. Another procedure for solving 
the problem of 1/s divergences was proposed in J18[ in the context of the U — oo out -of-equilibrium Anderson model: 



its method was to choose the zeroth-order spin Keldysh function (see section [II B 5 ) in order precisel y to cancel the 
divergences at order J 2 . This is completely equivalent to our approach, as explained in section [II C 2. 



Finally, we note that one may give a simple semiclassical derivation of the out-of-equilibrium result (0) based on 
a master equation. This derivation is due to L. Glazman and A. Kaminski Jl9| , and is presented in Appendix |A|. 
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The success of such a semiclassical approach (at this lowest order) is related to the fact that one can compute the 
equilibrium magnetization of a free quantum spin using a classical Ising model. In fact, at this order, the semi- 
classical master equation is strictly equivalent to the Keldysh component of the steady-state Dyson equation, so the 
two apparently disparate derivations yield the same result. 



B. Technical preliminaries: the Keldysh method 

Let us now turn to the technical details. In the following, for simplicity, we will write some equations for a generic 
fermionic held ip, which will be specialized afterwards to the fields representing the electrons and the spin. 



1. Generalities 



The basic idea of the non-equilibrium Keldysh method |2(J, |2l], £2| consists in taking the system at an initial 
time t = in an initial state described by a density matrix p and letting the system relax, using the Hamiltonian 
evolution given by H , to a long-time regime. In order to ensure that the system relaxes, it may be necessary to add 
some additional coupling terms to the Hamiltonian, in particular to break conservation laws: see the discussion in 
Appendix |e| Depending on the system, the long-time regime can be an equilibrium state, a non-equilibrium steady 
state, a non-time-translation-invariant steady state, or even an aging regime (in glassy systems). In the quantum dot 
problem, we assume that the system reaches at hnite voltage V a non-equilibrium steady state, in which we want to 
compute physical quantities. 

A "Keldysh" average of any quantity A is defined by 

{A(t)) K = (e iHt Ae- iHt ) Q = Tr (p e lHt Ae~ lHt ) , (11) 
where (. . .) is the average taken using the initial density matrix of the system. The steady-state average is given by: 

(A)=l\ra(A(t))^. 

Using the usual representation of the evolution operator e~ tHt as a T-ordered exponential in the interaction picture 
(and the anti-T-ordered one for e lHt ), one can obtain an expansion in the coupling constant J. Following the usual 
conventions, it is convenient to keep track of the two exponentials using a closed time contour, running from to +oo 



and back to 21 : we denote by + the upper contour (from to oo), which arises from expanding e lHt , and by — 
the lower contour. In accordance with this notation, we define the four Green's functions: 

G++(t,t') = -i(T1>{t)tf(t)) K , G-+(t,t>) ee -i{i,{t)^{t')) K , 

G^~(t,lf) = -i(Tfl>(t)^^)) K , G+-(t,t')= i(^{f)m) K - (12) 

Here T is the time-ordering operator, T the anti-time-ordering operator, and ip is any fermionic field. The two indices 
of the matrix G will be called "indices in Keldysh space". In the following, the equations (|l2| ) will be summarized 
with the notation: 

G^{t,t') = -i{{^{t)^{t'))) K (13) 
and G will always denote a 2 x 2 Keldysh matrix. 



The Green's functions (12) are not independent but can be expressed as functions of the retarded, advanced and 
Keldysh Green's functions defined respectively by p9| : 

R^t,t l ) = -iO{t-t')({^{t)^\t')}) K , 
A4t,t') = i8(t' -t)({^{t),^{t')}) K , 

K^{t,t') = -i([^{t),^{t')]) K , (14) 
as is shown from the transformation p2fl : 

c * s {%%~-) Mt^) MO-iM-M)- (15) 
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We will denote the first set of Green's functions ([L2|) "the ± basis" , and the second fll4| ) the "Larkin-Ovchinnikov 
(LO) basis". 

There are further relations between these Green's functions. In general, R^,{t,t') = (A^(t' ,t))* . Moreover, in 
equilibrium FDT reads: 

K^(oj) = h cq (aj)(^A 4 ,(uj) - #v>M) h eq (oj) = -tanh (j^) (16) 

thus the retarded Green's function is the only remaining independent Green's function, and so contains all the 
information about the state of the system. Out of equilibrium, however, K^,{oj) and should be independently 

determined. 

In order to write down the diagrammatic expansion, we also define the "bare" counterparts of the Green's functions 

in 



dig y) by: 



where the field V'int is in the interaction picture. From the expansions of the evolution operators, we obtain the 
diagrammatic expansion, provided that (. . .) satisfies Wick's theorem. For the problem at hand, we take as the 
initial condition the density matrix of the model with Jr = Jl = Jrl = and for the spin a finite s (arising from the 



coupling of the spin to a thermal bath; see section [II A above). Going to the long time limit, every function becomes 



a function of the difference of the times, and transforming to Fourier space we obtain standard Feynman rules. 

2. Lead electrons 

Let us begin with the lead electrons. We care here only about the local Green's functions, so we will drop the 
spatial indices. We denote by (Gcja^^^f,^) the Green's function describing the creation in lead /3 of an electron 
with spin er' at time t' ', and a corresponding annihilation in lead a of an electron with spin a at time t: 

{G c ) aai p a ,{t,t') = -i((c aa (t)cl,(t'))) . (17) 



(Here, the indices a, (3 £ {L, i?}, while a, a' = +, — , and the local electron operators are defined by c aa oc ^ k c^ko--) 
The Green's function G c is thus a priori a 4 x 4 matrix (in lead-spin space) , whose entries are themselves 2x2 Keldysh 
matrices. In the following, a bold notation with subscript c always designates such a 4 x 4 matrix. The bare Green's 
function G c o is however diagonal. The bare (diagonal) density of states is : 

{pco)ao,f3a' = + (tB c — II a )5 a /35 aa i 

(the energy levels are the same is both leads, but shifted by the Zeeman energy and the voltage). We assume that 
the baths, being much bigger than the impurity, are permanently in thermal equilibrium so that the bare electrons' 
functions read: 

RcoM = / de Pc ° { "\ + , (18a) 
J uj — e + iQ + 

K c o(^) = 2i7rh c0 (w)p c o(w), (18b) 

(h c o)acr,/3cr' = ^cq (^ — fl a ) S a p5 aa i , (18c) 

where the first two equations are matricial, and /i Q is the potential of the lead a — R, L. The voltage difference 
between the leads is given by V = /j-l — PR- 

3. Spin 

Since the spin operator is not appropriate for diagrammatic computations (it does not satisfy Wick's theorem) , we 
represent the spin 1/2 by three Majorana fermions r) a , a £ {x, y, z} which satisfy Wick's theorem and the relations 

S a = ~ l -e abc ifv c , (19a) 

fa')* = ?/\ (19b) 
{r, a ,r, b } = S ab . (19c) 
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Using (p9[), one can easily show that S satisfies the correct commutation relations and that S 2 — 3/4. Note that 
this last constraint is automatically satisfied, unlike in the case of a Dirac fermion representation, where a Lagrange 
multiplier would have been required to fix the magnitude of the spin: the Majorana representation therefore makes 
the computation simpler. In this paper, we consider only spin 1/2, but our computations could be extended to higher 
spin, provided that one used another representation for S. 

Let us now discuss the propagators G a b (a, b £ {x, y, z}) of the 77, in the presence of a magnetic field B s along the 
z direction. The general form is: 



G = I Gy X Gyy , G ab (t, t' ) = -% {{^(t)^^))) , (20) 




where the elements are 2x2 Keldysh matrices. Indeed G xz = G yz — by symmetry. To prove this, note that the 
Hamiltonian is invariant under a 7r rotation around the z axis, which is implemented by (j] x , rj y , rf ) — > (—r) x , — ?y y , rj z ) 
(and the corresponding rotation for the c electrons). Furthermore, making a n/2 rotation around the z axis, imple- 
mented by {rf , rj y \rj z ) — > (t] v , —rf,rf \ we obtain the relations 

G xy — Gyx: G X x — Gyy- (21) 

The Hamiltonian is also invariant under a 7r rotation around the x axis {{rj x ,rj y ,7y 2 ) — > {—T] x 1 r) y ,ri z )) together with 
a change of sign of the magnetic fields B s and i? c , and hence G xy and G^^ are respectively odd and even in the 
magnetic field. In particular, for B s = B c = 0, G xy — 0. The foregoing arguments apply to the full propagator and 
to the free propagator G° (computed with only the magnetic field B s ). 

A different basis is also useful: defining the Dirac fermion / and its Green's function Gb by 

/ - ^f, (22) 

G B (t) = - l ((f(t)p(0))) K , (23) 

we have the relations: 

G xx (lj)= G vv (u) = ^(Gb(w)+G-b(u>)), (24a) 

G X y{uj) = -G vx (uj) = % - (G b (lu) - G_ B H) . (24b) 

In the (/, rj z ) basis, the propagator is diagonal, so it is more convenient for example to write the Dyson equation, 
whereas the original basis {rj x ,r] y ,r] z ) is more convenient for the diagrammatics. The bare propagators in the (f,r] z ) 
basis are given by: 

R» = — (25a) 
u> + ts 

K° z (lu) = 2ift oq ( w )- 5 i- 5l (25b) 

R%(uj) = (25c) 

U) — B + is 

K B {u) = 24 q M (M _^ +a2 , (25d) 



where s is a small regulator which, as discussed in section |III A| above, should be thought of as the width due to 
coupling to an auxiliary thermal bath. 



4- Vertex factors and Dyson equations 

The vertex factors can be extracted simply from the Hamiltonian, and the Feynman rules are summarized in Fig. |l| 
(in the ± basis). Note that we have oriented the Majorana fermion lines, despite the operator property that = rj\ 
a proof that the lines are orientable is given in Appendix 

We now derive the expression for the Dyson equation describing the rj fields. In general a Dyson equation reads 

G-^lo) = G,T 1 H-£H, (26) 
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(G c o) QM , QM (6(j) 



T]2 _ T)l 



J, 




FIG. 1: Feynman rules: we represent the Majorana fermions with dashed lines, the electrons with solid lines, a, j3 are lead 
indices, u, v are spin indices for the electrons, k is a Keldysh index k = ±, rji, a € {x, y, z} are flavor indices for the Majorana 
fennion, and a are the Pauli matrices. The vertex is given in the ± basis. 



where the inversion has to be taken in the tensor product of the x, y, z space and the Keldysh space (thus with 6x6 
matrices). The free propagator is given by (p5|). In the x, y, z basis, the self energy can be written as: 



EH 



E d (w) iE (w) 
-iS a (u) S d (w) 
E,(w). 



E«(w) £*(w)\ , , , 



Using the definitions fl24|), we find: 



G 



±B 



I 



Eg B H + 2ish c(1 (m) \ 



T B - E^ B (ui) + is \u B - E% B (u)) + is\ 

1 







TB-£^ B (w)-is / 



(27) 



E± B (w) = E d (w)±E„(w). (28) 



C. Perturbation expansion in the steady state 



After these preliminaries, let us now come back to the perturbative expansion in the steady state and develop on 
the points introd uced in III A. In paragraph III C 1 we show that the limits J — > and s — ► do not commute; in 
paragraph III C 2 we solve that problem and present the derivation of Eq. (^) together with an algorithm for computing 
higher order terms in the perturbative expansion. 



1. Non-commutativity of limits 



We are interested in the steady-state values of the magnetizations (|2a|) and (|2b[), which can be expressed in the 
following way: 

„• r a, , / 

(29a) 
(29b) 



Pleads = - f — (— l) U (K c ) auou , 

Z I Ztt * — ' 



u = + , 



Indeed, the dot magnetization is given by: 



M dot = (5«) = -| (fo*,^]) = ^(t = 0). 
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-Micads is derived analogously. 



We begin with the expression for the dot magnetization derived from (29a) using the Dyson equation (|28| 



A/dot = 



4tt 



du> Im 



h eq (uj) 



u - B - + is 



dto £§ (w) - h cq (u)[^(u) - Eg (u/)] 



8i?r 



(30) 



where the last term denotes an antisymmetrization in B. The first term in ( p0|) can be transformed into a sum over 
Matsubara frequencies and thus we see that the limits J — > and s — > commute in that term. However, the second 
term is more interesting. Firstly, it vanishes in equilibrium since the numerator cancels, as required by the FDT. 
Moreover, using (for e — > 0) 



1 



B 



U{u>-B), 



(31) 



we see that in that second term these limits do not commute: 



• J — > and then s — > 0. We first take E^ = in the denominator. We find that this second term in (|3£ 
gives a diverging term at second order in J proportional to J\ L j 1 s. More precisely, using the expression for the 
self-energy at second order derived in Appendix ||, and ( |3l| ) with e = s, we find 



M, 



dot 



^tanh 



B 
2T 




up to a finite (i.e. not diverging as s — > 0) term of order J 2 and to 0(J 3 ) terms (ip is defined in (0)). 
• s — > and then J — * 0. Using (|3l] ) with e = ImE§(B) (we work at dominant order in J), we find: 



A/dot = — 



1 



4 VE£(£)-E£(£) 



- (B -B) + O(J) 



(32) 



where the self-energy has to be expanded at order 2 in J (the first term in (|3^) cancels a part of the second 
term). 



The second limit is the physical one (as explained in section II) and it gives a formula (32) for the magnetization at 
order 0. However, to make use of this formula one needs to know the Keldysh (£f (B)) and spectral (E|(B) - £§(£)) 
parts of the self-energy. The leading terms of these self-energies are of order 2 in J, but because the ratio of them 
is taken, they determine the magnetization at order 0. The crucial point is that d32| ) is in fact an implicit equation 
for hs{B), which appears on the left-hand side because Afdot is defined in terms of it by (|29|), and on the right-hand 
side because the self-energies at 0(J 2 ) depend on the Keldysh Green's function at O(l), and hence on hs- The most 
elegant way to capture this 'feedback' effect is to reformulate the problem in terms of the perturbative expansion of the 
steady-state Green's function in powers of J, starting from the Dyson equation in which the s — > limit has already 
been taken. We now describe this method in more detail, and specify the procedure for calculating the steady-state 
Green's functions perturbatively in J to arbitrary order. 



2. Perturbative expansion of steady-state quantities 



Let us now give a general method for computing the perturbative expansion of the Green's functions, and use it 
to derive Eq. (0). First, we reformulate slightly the diagrammatic expansion in terms of the full Green's function 
and of skeleton diagrams, and we explicitly derive the steady-state Dyson equation as a functional equation for 
(again, for simplicity, we write some equations for a generic field ip). The "generic" Dyson equation reads: 



G 



{J,s,u) = (g^(uj) - £ skc i[Gy(J,s,w) 



(33) 



where E s k e i is a functional of G^ defined by the skeleton expansion of the self-energy (see Appendix ^ for a derivation) . 
We first take the s — > limit in (33) and then solve it order by order in J. It is in principle sufficient to solve (|33j) 
for Gi/, since the magnetization can be directly extracted from a Green's function, and other physical quantities (e.g. 
currents) are given by their skeleton expansions. 
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We first derive the explicit form of the Dyson equation for all fields. Defining 



K 



yK 



yR 



V 2VK 2l7T 

the full set of Dyson equations for the electrons and the Majorana fields can be rewritten explicitly as: 
TL c (u>) 



K c (w 

Rb(u 

h B (uj 
R z (lj 
K z {uo 
h z (u> 



1 



R. 



(l-R c oSfi) _1 K c0 (l 
1 

w- J B-£g(w)+iO+' 
h B {uj)\A B {uj) - R B {u) 



_ U B 



a B (oj) 



1 



w-S^(w)+iO+' 
h z (u>)(A z (u)-R z (w] 



a z (uo) 



(34) 

(35) 

(36) 
(37) 

(38) 

(39) 

(40) 

(41) 

(42) 



where the bold symbols are 4x4 matrices (in lead-spin space). For completeness, we have also written the definitions 
of the h functions. We have three blocks of equations, for the electrons, the / field and the rf field respectively. 
Within these blocks, we have an equation for the retarded function ( |35| , ^7| , p0 ) , a n equation for the Keldysh function 
( |3l|j^fill ) , and for the Majorana fermions the definition of the h function^ 39| , [42"|) . We define h B with ( |39|) rather 
than with (|38|) since the spectral density is a delta peak at order 0, whereas the self-energy is a smooth function. 

The spin and the lead electrons appear on a different footing: the order part of the electronic Keldysh function 
is given by K c o whereas the order parts of h B and h z must be computed using equations ( ^Qp^ ). Note that to 
obtain the order n part of h B , one needs to compute the self-energies at order n + 2. The method of obtaining the 
perturbative expansion in J order by order is as follows (denoting by the order n part of any function /): 

1. Assume that we have the expansion of all functions to order n — 1. 

2. Since S s k c i is at least of order 1, using (^5|,^7]|4C|) , we compute R^ B \ i?i"' > , 



R c (n) and K c (n) . 

3. Since S^ c | and Sg kc i begin at order 2, we compute the a K and a to order n + 2, as functions of the unknowns 



^skcl 

hP and h z n) . 

4. We then obtain closed equations for hg and h z from 



3) (expanded to order n). 



Thus the order part of the impurity magnetization is given by (from (29a)): 



Mdot = 



(B). 



The order parts of R B and K B are 



R 



(0) 



1 



K ( ° } = 2iirh { ° ) (B)S(lu - B), 



(43) 



(44) 



and the bare Green's functions of the electrons are given by ( Jig} . We compute h B (u> — B) at order by expanding 
the self-energies at second order (given by the diagram of FigTg), and then solving (|39|) for h B (B). Finally we find 
(Q) in the large bandwidth limit (D — > oo). The computation is presented in detail in Appendix [f]. This completes 
the computation of the Green's function to order 0. 

Note that K c is thermal at order 0, which expresses the fact that the leads are in thermal equilibrium. Therefore 
the leads' contribution to the total magnetization at order is given by the Pauli term, which was explicitly excluded 
from A/tot! hence M to t = M^ot at this order, as claimed in (13). Moreover, we have not written explicitly the full 
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a 




FIG. 2: Diagram of the skeleton self-energy at second order. 

forms of the functions and h z (uj) to this order, since they are not required in the calculation of the zeroth-order 

Keldysh functions: Hb(B) and h z (0) = are sufficient, since the spectral density (at this order) is a delta peak. The 
full functions would, however, be needed to compute at second order; the function /ib(w) can easily be extracted from 
Appendix |], and the calculation of h z (u>) proceeds along similar lines. 

Our earlier interpretation of the l/s divergences in 'straightforward' perturbation theory is borne out by this result. 
As stated above, the l/s divergences result from an incorrect choice of zeroth-order distribution function h. If we 
insert the 0(1) part of ( |39|) into (^o|) in place of h cq , we see that the divergences are cancelled at order J 2 , since to 
this order we are now using the correct long-time distribution function. This shows that our method and that of |I§| ] 
are equivalent. 

It is important to note that these corrections to the zeroth-order terms in perturbation series are in no way restricted 
to quantities such as the magnetization. On the contrary, since what we have really calculated is the correction to the 
zeroth-order Keldysh Green's function, they manifest themselves in many quantities. As an example, we may consider 
the current-current correlator: the leading terms in this quantity are of order J 2 , and are calculated by inserting the 
zeroth-order Green's functions into the skeleton diagrams shown schematically in Fig. 0. The result is that 




FIG. 3: The two skeleton diagrams contributing to the current-current correlator at leading order. In each diagram, one of 
the electrons (represented by the solid lines) is from the left lead, and the other is from the right. The vertices represent I, the 
current operator. 



(I(-oj)I(oj)) = 



7T {PqJlrY 

2h 2 



hf (B)] p(w, V) - hf (B) [p(co + B,V)- p(u - B, V)] 



+ \p(cj, V) + P (uj - B, V) + p{lo + B, V)] , 



(45) 



where the function p(x, y) is defined by 



p{x,y) 



x-y 



x + y 



aP(x-y) _ I e f3(x+y) _ \ ' 



(46) 



The result ( |45| ) is a function of hg)(B), and is therefore clearly sensitive to the corrections made to the zeroth-order 
Green's function, as expected on the basis of the discussion above. 
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IV. CONCLUSION 



The purpose of this note has been to present the expression for the order magnetization in the Kondo model out 
of equilibrium (Eq. Q), and a systematic procedure for obtaining higher-order corrections to this result. The result 
we obtain may seem surprising, in the sense that even at order in J it does not coincide with the equilibrium 
expression | tanh(£>/2T). Indeed the out-of-equilibrium distribution function which describes the population of the 
two levels of a weakly coupled spin is in general not thermal but must be computed by solving a transport equation: 
it is determined by the steady state into which the voltage difference forces the system. Moreover, this distribution 
function also enters the computation of other physical quantities (e.g. currents) and their perturbative expansions 
therefore exhibit similar phenomena. 

Finally, we stress that the issue discussed in this note is not directly related to the so-called "decoherence time" 
issue. Answering the question of strong coupling at finite voltage requires computations at higher orders in J. It is 
clear from the above, however, that the behavior of the out-of-equilibrium perturbative expansion will be markedly 
different from that of its equilibrium counterpart. For example, the O(J) contribution to the Keldysh Green's function 
of the Majorana fermions depends on the 0(J 3 ) contribution to the self-energies. These, however, contain terms that 
diverge like lnD in the D — > oo limit, and so such logarithmic divergences may be expected to appear at O(J) in some 
of the Green's functions, and therefore in physical properties such as the magnetization. Indeed, this phenomenon has 



been reported recently 1 23 ; the full interpretation of this striking departure from equilibrium behavior merits further 
work. 
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APPENDIX A: A SEMI-CLASSICAL DERIVATION OF (@) 

The magnetization at order in J (0) can also be derived using a semiclassical master equation approach |l^| . By 
symmetry, the reduced density matrix of the spin is diagonal, so we consider the spin as classical and characterize its 
state by the probability of its being up, P-j-, or down, P^. The spin dynamics is governed by a master equation: 



dt 
dP L (t) 
dt 



riTnM-r u nW, (Ala) 
r n p T (t) - r iT p ; (t), (Aib) 



where Tab is the rate of the spin-flip process A — > B induced by the Kondo terms. At second order in perturbation 
theory, these rates are given by: 



r lT ~ J "/3 / *i de 2 ^ B c - V>a)p( e 2 + B c - V>p)riF{e\ - McOU - n F^2 - - e 2 + B s ) (A2a) 

r Ti ~ ^2 J a/3 J deide 2 p(ei + B c - p a )p{e 2 - B c - (ip)n F (ei - p a )(l - n F (e 2 - n l3 )^S(e 1 - e 2 - B s ) (A2b) 
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Taking the large bandwidth (D — > oo) limit and integrating, the rates become: 

1 Tl - l/i? + J Lj e B s /T _ 1 +J RL\ e ( Bs -V)/T _ 1 + e (B s + V)/T _ \ I ' \ Aa& ) 



^ii + Jh) e - Bs /T _l +J RL\ e -(B s +V)/T _ l + e _(B s -y)/T _ 1 J ' [ 6 ' 



The steady state value of the probability is given by equating the left hand sides of (Al) to 0, and the magnetization 
is given by M dot = (Pf - P[)/2, thus 

M d ot=i(&^>y (A4) 



2 



IT f 1 Tl 



which leads to Eq. (||). 

One can see that the Dyson equation in steady state, at second order in J, maps exactly to the master equation, 
although they appear to have different transients. So at order 0, this computation is just a reformulation of the one 
presented above. 



APPENDIX B: A FREE MODEL 



In this appendix, we recall [|14| |2J|] the solution of a simple free model in an out-of-equilibrium steady state regime, 
which displays a result very similar to (^). We consider a free level coupled to some reservoirs of free electrons via 
hopping terms. The Hamiltonian is given by: 



H = ^ £kC^ k c Qk + ^ 9a (eld + d)c a 



(Bl) 



where the local reservoir-electron operator c a is defined by c a — c a (x = 0) = SI -1 / 2 ^ k c a k- (fi is the usual 
normalization factor related to the volume of the reservoirs.) As before, we assume that the reservoirs are thermalized 
with the same temperature as each other, but with different chemical potentials /i Q . In particular, we are interested 
in a model with three reservoirs: 1 and 2 are the leads (at different potentials) and 3 an additional thermal bath (with 
fis = 0) to which the level is coupled. Since it is a Gaussian model, we can simply solve the Dyson equation: 



G- 1 = Gq 1 - E. 

We use the (d, c\, c%, C3) basis. The inverse bare Green's function is given by: 

G 1 H = 



(B2) 





- £d 








° ^ 




























{ 














(B3) 



where G a is the bare Green's function of the reservoir a, given in the large bandwidth limit by (in the LO basis): 



G-\lo) = — 
ttPo 



i -2ih cq (u> - Ha) 
-i 



(B4) 



where po is the density of states, which we take to be the same for each reservoir. The self-energy is given by: 



EH 



/0 91 g2 53 N 

9i 

52 

V53 0, 



Solving (jB2|), we find the occupation of the dot in the steady state: 

n(e d ) = (dU) . 



(B5) 



(B6) 
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In the limit where 31,52,53 go to zero in fixed ratios, we obtain: 

/ \ 9in F (e d - fix) +gln F (e d - /i 2 ) + 5 |n F (e d ) 

n{ed) = gj+gl + g! (B7) 

where n F is the Fermi function. The properties of this result are similar to those of (^): the occupation of the dot, in 
the limit of zero couplings, depends on the ratios of these couplings and is not in general given by the Fermi function. 
Moreover, if 33 3> g%, 32, we recover the equilibrium result since the physics is dominated by the thermal bath 3. If, on 
the other hand, we take 33 31 , <?2 , we find a non-equilibrium result since the occupation at order in perturbation 
theory in the g's is determined by the leads 1 and 2. 



APPENDIX C: ORIENTABILITY OF MAJORANA FERMION LINES 



In this appendix, we demonstrate that we can treat the Majorana lines as oriented in the diagrams. To show this, 
it is simplest to take the lines to represent the bare zero-field Majorana functions, and to treat both the magnetic 
field and the Kondo interaction as vertices. We have three species of Majorana fermion, {rj x ,rf ,rj z }; formally, we 
may represent each of these as the sum of the creation and annihilation operators of a Dirac fermion: 

v a = ^ (/. + /*); (ci) 

the dual operator f a — /J decouples from the problem and may be ignored. In the absence of a magnetic field and 
interactions, all Majorana Green's functions are diagonal, so the only Green's functions that occur are: 

G a = -i((v a (t)v a (0)))K 

'—((faWm))^ +p.h., (C2) 

where the notation "p.h." stands for "particle- hole" , i.e. / <-> p . We thus see that the bare Majorana Green's 
function may be written simply as the sum of two bare Dirac fermion Green's functions of opposite orientations. 

It is easy to see that the vertex factors at the magnetic field and Kondo vertices do not depend on whether the 
Dirac fermion lines are incoming or outgoing. In the case of the magnetic field, for example, the interaction term is 

= iBrfrf 

= Y (fx + ft) (fy + , (C3) 

so we see that the vertex factor is the same irrespective of the orientation of the two /-lines; a parallel argument may 
be given for the Kondo interaction vertex. Hence, each diagram consists of a sum of 2 N diagrams that differ only in 
the orientations of their N Dirac fermion lines. 

But these orientational differences do not alter the value of the diagram, since the bare Green's function of the 



/-fermion is particle- hole symmetric, and hence (C2) corresponds to the Green's function for a single orientation of 
the Dirac fermion line, with the prefactor 1/2 removed. Hence we may represent the Majorana Green's functions in 
all diagrams using oriented lines. 



APPENDIX D: ERRATUM TO juj 

In a previous paper |l5| , a formula (Eq. 2) was proposed for the (second order) perturbative expansion of the 
magnetic susceptibility in the out-of-equilibrium steady state which was of the form: 

X(T,V) = ^ + 0(J). (Dl) 

In this appendix, we briefly rediscuss its derivation and explain why it is incorrect. The method used was straightfor- 
ward Keldysh perturbation theory to second order with s finite in which the bare Keldysh function of the Majorana 
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spin was taken to be thermal. The l/s divergences were regulated using a "point splitting" procedure: since they 
occur due to the coincidence of two poles in the integrals, one splits these poles on the real axis to a distance 8 and 
sends 5 to at the end of the computation, term by term. The justification given in footnote 20 of |ljj is however 
incorrect: endowing the Majorana fermions on the dot with a fictitious dispersion does not lead to this prescription, 
since divergences reappear when the bandwidth is sent to zero. The "point splitting" rcgularization prescription can 
however be described physically as follows: let the applied magnetic fields B s and B c oscillate slowly at a frequency 
ujq; compute the total magnetization at frequency luq up to second order in J; then take the limit loq — > in the 
coefficient of each power of J. 

However, this regulation procedure is based on an interchange of the order of the limits ujq — > and J — > 0. We 
want to calculate the static magnetic response, and consequently wish to take lo — ► before J — > 0; but in fact the 
technique used in |Iq ] does the opposite: it expands in J ( J — > 0) before taking the ljq — ► limit. It is simple to show 
that these limits do not commute; this can be seen explicitly from the form of the second term of ( |30| ) : 

Re / dL0 < i r^rr rv (D2) 

where g is some function and 7 ~ £"(£?) (the limit s — > having already been taken). If we take the limits J — > 
(i.e. 7 — ► 0) and then wo -> we obtain the result of the point splitting prescription of fl5| : 

du ( 9 -M] , (D3 ) 



whereas if we take them in the opposite order, ujq — > and then J — > 0, we recover (^). 



APPENDIX E: THE STEADY STATE DYSON EQUATION 

In this appendix, we present a formal derivation of the steady state Dyson equation ([33]). First of all, let us 
emphasise that our calculation is carried out in the steady state regime, that is in the long-time limit after the 
switching on of the interactions. We assume that, in this limit, the system evolves to a time- independent steady 
state under the time evolution described by its Hamiltonian (as required by the Keldysh method: see section III B|) . 

i 



Strictly speaking, this is not possible, since the Hamiltonian (hi) conserves the total magnetization of the system and 
that conservation law prevents the magnetization of the system from relaxing. However, this conservation law is not 
physical (we have omitted, for example, spin-orbit terms in the leads); therefore, to allow the system to relax to its 
steady state, we introduce a coupling g that breaks the conservation laws. As a specific example, one could consider 
an anisotropic (J x =/= J y =/= J z ) Kondo model. As this extra coupling is relaxed to zero, the transient time taken to 
reach the steady state diverges but we make the assumption that the values of physical quantities in the steady state 
are smooth functions of g. Therefore, once we have taken the t — > 00 limit, we can set g = in the equation which 
determines the steady state. 

The details of the derivation are as follows: 

1. Using the perturbative expansion, we establish the Dyson equation at finite times: 

<lu (G^(t,u) - E ske i[Gy (t, u)j * G 4 ,{u,t') = S(t - t') ® 1, (El) 

where S s kei is a functional of defined by the skeleton expansion of the self-energy. The product should be 
understood as a matrix product in the LO basis. Gj> is a function of two times and of J, g. To obtain this 



equation, we write the Dyson equation in the finite time dia gram matic expansion outlined in section III B 4 j22| 



and use the definition of the skeleton diagrams [|26| , |27| . In (El), the times t, t' and u run from —00 to 00, and 



the couplings J are time dependent: J a (t) = J a Q(t), i.e. we switch on the interaction suddenly at t = 0. 

2. The assumption that the system relaxes to a non-equilibrium time translation invariant steady state (the coupling 
g to the relaxation bath is finite) is transcribed mathematically as the existence of the limit 

lim G(t,t')=G(r) (E2) 



(we shall denote it with the same function name). This assumption is not trivial. In particular, in |28|, the 
existence of steady oscillating states has been suggested using a large- N slave boson treatment. We exclude them 
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here on physical grounds, since we do not expect such states to appear in the regime where the perturbation 
theory is applicable anyway (i.e. at high temperatures or at high magnetic fields in the Kondo problem). We 
can thus take the long time limit, and Fourier transform the Green's functions and self-energy in flEl| ) to obtain: 

G^(J,s,g,u) = (g^q(cj) -T, ske i[G^](J,s,g,u)^ . (E3) 

3. Taking the g — ► limit as discussed above, we finally obtain the Dyson equation (p3|): 

G^{J,s,u) = (g^q(u) - S ske i[Gv,](J,s,w)) _1 . (E4) 

APPENDIX F: COMPUTATION OF THE SELF-ENERGY DIAGRAM 

In this appendix, we present the computation of hs(B) at order 0. The computation is in three steps. First, we 
compute the self-energy diagrams at second order as a function of a "generic" diagram. (This simplifies the problem 
by separating the Majorana and spin indices from the Keldysh structure.) Second, we compute this generic diagram. 
Finally, we solve the resulting implicit equation for hsiB). Throughout this appendix, we streamline our notation 
by omitting the temperature T; in the final formulas, therefore, B and V should be replaced by B/T and V/T 
respectively. 

The first part of the computation reduces the spin and lead indices, and thus expresses the self energy diagrams as 
functions of the "generic" diagram presented in Fig. |], where we allow any potentials V\ and V 2 for the electrons and 
any field B for the internal Majorana line. The main formula is: 




FIG. 4: Diagram of the "generic" self-energy at second order. 



E 



\ J ap\ 



<T=+,- 



f(u, B s , V a ,o-, Vp,a) + /(w, 0, V a , 1, Vp, -1) 



(Fl) 



In this expression, the Keldysh structure is implicit and by convention, the Majorana line is an rf line when B s = 
and an / line otherwise. To establish (Fl), we compute the spin and lead indices of Yi xx and Y, xy which are given by 
the diagrams of Fig. || (with the Feynman rules given in section III B ) , and we use (28). 
We now compute the Keldysh structure of the generic diagram /. We have: 



f+-(u) = m J dxdyp(x + a t B c - Vx)p{y + a 2 B c - V 2 )n F (x - Vx)h F {y - V 2 )(l + h{B))8{u +y-x-B), 
f-+(u) = -iir J dxdyp(x + a 1 B c -Vi)p(y + a 2 B c -V 2 )n F (x-V l )nF(y-V 2 )(l-h(B))5(uj + y-x-B), (F2) 
where n F (x) = n F (—x), and the function h(u) is defined by 

h{uo) 



We then use the relations |22| 



h B {oj) w^O, 
h z (0) u = 0. 



f K = -(/+- +/-+), 
f A -f R = -(/+- -/-+)■ 



(F3) 



(F4) 
(F5) 
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FIG. 5: Diagrams ofS xx and ~S xy 
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FIG. 6: Computation of the Keldysh structure of the generic diagram. 



We can perform the integrals in the large bandwidth limit D — > oo using 

A 



dxn,F{x + A)np(x) 



e A -l 



and replacing the densities p by their finite value po (we can take the limit D — > oo under the integral). Using the 
definition V = V2 — Vi, we find: 



(f A -n 



2iir 



1 



tanh 



u-B + V 



KB) 



(F6) 
(F7) 



Using (Fl), and introducing 



(p(x) 



tanh § ' 



(F8) 



we now have: 



\PoJrl\' 



u b {oj) = , ' h B {B)(uj - B + V) - h z (0)(u> + V) + ip{u) - B + V) + <p(uj + V) + (V -> -V") 
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of( W )=- 



E 

a=R,L 

\poJrl\ 2 



\PoJa 



2 r 



-h B {B){uj - B) - h z (0)(oj) + <p(u -B) + v (u) 
h B (B) 



(F9) 



(u -B + V) [1- 



tanh^|±^ 



E 

a=R,L 



\poJa 



{u-B)[l 



h B {B) 

j-B 



+ (w + V) ^1 

ft*(0) 



ftz(O) 



tanh^ 



+ (U^-V) 



tanh 



+ w 1 



tanh | 



Using the definition of the Majorana Green's function G^ b (t), 

G% > (t) = -i{[r h (t)M0)]) 



(F10) 



(Fll) 



we have G^ b (—t) = —G ba {t). In the frequency representation, this reads G^ b (uj) = — G ba {— w) and hence we infer that 
Gf z (uj) is odd in frequency. Consequently, given the form G^ z — h z (ui)8(oj) it is clear that 

h z (0)=0. 

Using this result together with ( |39| ) and denoting x = —hs{B), we have: 

\Jrl I 2 (2x<p(V) + 2B) + (4 + J| ) (2x + B) 



\Jrl\ 2 (2<p(Y) + <P(B + V) + y{B - V)) + (J% + J|) (2 + <p(B)) ' 
Solving for x and substituting into gives Eq. (|) of the text. 
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